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ABSTRACT 


In  this  thesis  a  feasibility  study  is  presented  for 
the  measurement  of  thermal  and  normal  momentum  accommodation 
coefficients.  Two  molecular  number  density  detection  sys¬ 
tems  are  analysed:  an  ionization  detector,  and  an  electron 
flourescence  device. 

An  analysis  of  experimental  accuracy  was  carried  out 
on  the  expression  for  the  thermal  accommodation  coefficient. 
A  one  percent  standard  error  in  the  measurement  of  the 
parameter s  wa s  ass umed  to  determine  which  measurements  were 
most  critical.  The  measurements  which  were  found  to  be  most 
critical  were  the  number  density  in  chamber  A  and  the  number 
density  in  the  active  zone  of  the  detector. 

The  instrumentation  was  selected  to  minimize  the  error 
associated  with  the  measurement  of  the  most  critical  para¬ 
meters  . 

The  percent  most  probable  error  in  the  measurement  of 
the  thermal  accommodation  coefficient  is  less  than  six  per¬ 
cent.  The  percent  most  probable  error  in  the  measurement 
of  the  normal  momentum  accommodation  coefficient  is  less 
than  five  percent. 

Both  systems  are  well  suited  for  measuring  these 
accommodation  coefficients. 
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CHAPTER  I 


INTRODUCTION 

Free  molecular  flow  conditions  exist  when  the  mean 
free  path  of  molecules  in  a  gas  is  large  compared  to  a 
characteristic  dimension  on  a  solid  surface.  Free  molecular 
flow  conditions  exist  in  the  atmosphere  above  an  altitude 
of  approximately  80  miles.  In  this  ambient  atmosphere,  low 
orbiting  earth  satellites  experience  orbital  decay  due  to 
the  exchange  of  momentum  between  the  atmospheric  gas  mole¬ 
cules  and  the  satellite  surface. 

One  of  the  earliest  and  most  successful  techniques  for 
the  measurement  of  thermal  accommodation  coefficient,  a,  is 
the  conductivity  cell  method  first  introduced  by  Knudsen  in 
1910.  The  macroscopic  heat  transfer  by  free  molecular  con¬ 
duction  between  a  central  electrically  heated  filament  and 
a  surrounding  wall  is  measured  directly  to  obtain  the  accom¬ 
modation  coefficient  for  the  filament  material. 

Roberts  was  apparently  the  first  to  realize  how  import¬ 
ant  the  effects  of  surface  contamination  can  be,  and  in  fact 
all  measured  values  of  accommodation  coefficients  are  reli¬ 
able  only  to  the  degree  that  the  condition  of  the  solid  sur¬ 
face  is  known.  One  of  the  advantages  of  the  conductivity  ' 
cell  method  is  that  this  relatively  simple  apparatus  makes 
it  easy  to  obtain  clean  metal  surfaces  by  flash  heating  to 
evaporate  contaminant  gases.  Some  disadvantages  are  that  the 
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range  of  solid  surface  materials  is  limited  to  conducting 
metal  filaments,  and  that  the  solid  surface  temperature  range 
is  limited  by  radiation  heat  transfer  effects  which  become 
large  if  the  filament  temperature  is  much  above  the  tempera¬ 
ture  of  the  surroundings. 

The  conductivity  cell  method  and  results  are  discussed 
in  a  review  paper  by  H.Y.  Wachman  (1). 

More  recent  molecular  beam  techniques  used  by;  Datz, 
Moore  and  Taylor,  Smith  and  Fite,  and  others  are  described 
in  a  review  by  F.C.  Hurlbut  (2) .  This  method  is  concerned 
with  the  scattering  pattern  of  the  reflected  molecules.  This 
method  yields  some  microscopic  behavior  of  the  gas-surface 
interaction . 

In  this  thesis,  a  feasibility  study  is  made  on  a  tech¬ 
nique  for  measuring  thermal  and  normal  momentum  accommodation 
coefficients.  The  vacuum  system  used  for  the  experiment  is 
designed  to  control  impurity  gases  (thus  surface  contamina¬ 
tion),  and  is  designed  so  that  any  material  (metallic  or  non- 
metallic)  can  be  used  for  the  test  surface.  Radiation  effects 
do  not  limit  the  test  surface  temperature  range. 

The  main  concerns  of  this  thesis  are:  the  formulation 
of  a  complete  error  analysis,  design  of  a  device  which  moni¬ 
tors  number  density,  and  the  selection  of  instrumentation 
required  for  the  best  possible  accuracy.  Two  devices  for 
monitoring  the  number  density  will  be  considered;  an  ioniza¬ 
tion  detector  and  an  electron  flourescence  device.  The  two 
systems  will  be  compared,  and  the  required  system  adaptations 
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will  be  cited. 

1 . 1  A c c o mm odation  Coefficients 

The  normal  momentum  and  energy  transfer  at  a  gas-solid 
interface  are  the  quantities  which  can  be  determined  by  the 
experimental  apparatus  to  be  described  in  this  thesis. 

The  average  energy  flux  at  a  gas-solid  interface  is  the 
difference  between  the  energy  of  the  incident  molecules  and 
the  energy  of  the  reflected  molecules.  The  incident  mole¬ 
cular  stream  and  the  reflected  molecular  stream  may  be  con¬ 
sidered  independently  if  free  molecular  flow  conditions  exist. 
If  the  mean  energy  carried  to  the  surface  per  incident  mole¬ 
cule  is  (corresponding  to  temperature  T^),  and  if  the 

number  of  molecules  striking  a  unit  area  of  the  surface  per 
unit  time  is  Nj_;  then  the  total  amount  of  energy,  ,  carried 
to  the  solid  surface  is  £iNi .  Similarly,  the  amount  of  energy 
E r ,  leaving  the  solid  surface  is  £rNr  (corresponding  to  a 
temperature  Tr ,  where  T±  <  Tr  <  Tg  if  Ts  >  T±) .  Therefore, 
the  energy  flux  is  E^  -  E r .  If  the  gas-solid  interaction 

were  complete  the  temperature  of  the  reflected  molecules  would 
be  fully  accommodated  to  the  surface  temperature,  Tg.  Hence, 
the  energy  flux  in  the  completely  accommodated  case  is  E^  —  Eg 
These  two  energy  fluxes  suggest  that  the  average  effi¬ 
ciency  of  the  energy  exchanged  at  a  gas-solid  interface  can 
be  expressed  as  a  semi-emperical  coefficient.  This  coeffi¬ 
cient  was  defined  by  Knudsen  in  1910  as  the  energy  accommoda¬ 
tion  coefficient,  Ot. 
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(1.1.1) 


Similarly,  the  normal  momentum  accommodation  coeffi¬ 
cient,  0} ,  can  be  defined  as: 


c 


(1.1.2) 


Here  the  subscripts  have  the  same  meaning  as  in  expression 

1.1.1. 

The  above  expression  may  be  reduced  to  a  more  conveni¬ 
ent  form  by  assuming  the  velocity  distribution  of  the  reflec¬ 
ted  molecules  to  be  Maxwellian.  It  can  be  shown  that, 


(1.1.3) 


and  , 


a  ’  = 


( T  ^ ) 1 / 2  -  (Tr ) 1 / 2 

(T  .  ) 1 ^  2  -  (T  )1^2 
v  l  s 


(1.1.4) 


for  a  stationary  gas  in  thermal  equilibrium. 


1 . 2  Vacuum  System 


The  vacuum  system  used  for  the  measurement  of  a  and  o' 
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is  described  in  detail  in  reference  (3).  Modifications  will 
be  cited  in  a  later  chapter. 

For  the  sake  of  explanation,  a  schematic  of  the  system 
is  shown  in  Figure  1.3.1.  The  working  pressure  in  the  cham¬ 
ber  A  is  maintained  at  1  x  10~^  torr.  The  effective  pumping 
speed  from  chamber  B  is  1900  litres  per  second.  Chamber  A 
and  B  are  connected  by  a  2.08  mm  diameter  orifice  which  has 
a  molecular  conductance  of  0.37  literes  per  second.  From 
molecular  flow  theory,  the  pressure  in  chamber  B  should  be 
approximately  1.95  x  10“^  torr. 

The  instrumentation  will  be  cited  in  sections  2.5  and 
3.5  after  a  complete  error  analysis  has  been  carried  out. 

The  error  analysis  will  enable  us  to  determine  the  type  and 
accuracy  of  instrumentation  required. 

The  devices  used  for  monitoring  the  number  density  of 
the  molecular  beam  will  also  be  discussed  in  sections  2.1 
and  3.1.  , 

1 . 3  Test  Surface 

According  to  Wachman  (1),  early  investigators  failed 
to  specify  the  amount  and  type  of  contaminants  present  during 
the  test.  Thus  the  values  of  the  accommodation  coefficients 
reported  are  not  reliable.  To  overcome  this  problem,  a  method 
used  by  J.J.  Hinchen  and  E.  Sheperd  Malloy  (6)  was  employed 
to  clean  the  platinum  test  surface  prior  to  the  experimental 
run.  Polycrystalline  platinum  targets  are  resistively  heated 

This  procedure  cleans  the  surface  by 


in  ultra  high  vacuum. 
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oxide  decomposition  and  causes  the  growth  of  large  crystals 
with  the  < 1 1 0 >  plane  parallel  to  the  test  surface. 

A  platinum  -  10%  platinum,  90%  rhodium  thermocouple 
with  a  suitable  millivolt  readout  can  be  used  to  measure  the 
temperature  of  the  surface.  A  0.003  inch  wire  was  used  to 
avoid  large  heat  loss  by  conduction  along  the  thermocouple 
wire . 

The  temp erature  of  the  incident  gas  molecules  is  assumed 
to  be  equal  to  the  wall  temperature  of  chamber  A.  Chromel- 
alumel  thermocouples  together  with  a  suitable  potentiometer 
can  be  used  to  determine  the  wall  temperature.  The  wall  of 
chamber  A  can  be  kept  at  a  constant  uniform  temperature  by 
circulating  water  through  the  water  jacket. 

The  test  gas  chosen  was  nitrogen  because  it  does  not 
chemisorb  on  platinum.  Nitrogen  has  a  well  defined  (0,0) 
band  in  the  first  negative  system  and  is  therefore  well  suited 
for  the  electron  flourescence  technique. 
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CHAPTER  II 


IONIZATION  DETECTION 


2 . 1  Description 

The  number  density  of  the  molecular  beam  will  be  mea¬ 
sured  by  a  device  which  is  similar  to  the  detection  system 
designed  by  Hagena  e t.  a  1 .  (4)  at  the  University  of  Virginia. 
This  device  mainly  consists  of  three  parts:  (a)  ionization 
chamber;  (b)  ion  collector  and  (c)  filament  and  shield.  A 
schematic  diagram  of  this  detector  is  shown  in  Figure  2.1.1. 

The  ionization  chamber  was  machined  from  stainless 
steel.  The  aperture,  enabling  the  molecules  to  pass  through 
the  electron  film,  is  covered  by  a  90%  transparent  stainless 
steel  mesh.  This  mesh  provides  a  uniform  electric  field  to 
accelerate  the  positive  ions  towards  the  ion  collector. 

The  ion  collector  consists  of  a  stainless  steel  frame 
covered  with  90%  transparent  mesh.  The  collector  is  mounted 
on  the  ionization  chamber,  opposite  the  entrance  aperture. 
Teflon  is  used  as  the  insulating  material. 

The  filament  shield  arrangement  is  also  mounted  on  the 
ionization  chamber,  but  opposite  the  electron  beam  slits. 

This  arrangement  has  the  advantage  that  the  cathode  and  anode 
may  be  interchanged  by  external  switching.  Electrons  are 
produced  by  rss is tively  heating  the  0.010  inch  tungsten  fila 
ment.  The  filament  and  shield  are  insulated  from  the  ioni¬ 
zation  chamber  by  alumina  ceramic  beads. 
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The  electron  film  is  collimated  by  a  900  gauss  magnetic 
field.  The  permanent  magnet  supplying  the  field  was  enclosed 
in  a  stainless  steel  case  to  overcome  outgassing. 

2 . 2  Theor  y 

Since  the  orifice  connecting  chambers  A  and  B  of  the 
vacuum  system  has  a  diameter  which  is  much  less  than  the  mean 
free  path  for  nitrogen  at  a  pressure  of  10“^  torr  .  The  number 
density  at  the  detector  can  be  calculated  from  the  kinetic 
theory  of  gases. 

The  number  of  molecules  with  velocities  between  c  and 
C-  +  dc_  is  given  by  fdc^,  where  f  is  the  distribution  function. 
If  the  molecules  are  to  cross  an  area  dA  in  time  dt  ,  they  must 
have  occupied  a  volume  c^.n  dt  dA,  where  n  is  the  unit  vector 
normal  to  dA  (Figure  2.2,1).  Since  dc_  =  c ^  sin  <j)  d(f>  dG  dc, 
the  number  of  molecules  crossing  dA  in  time  dt  with  speeds 
between  c  and  c  +  dc  is  : 

c^  f  dc  d0  sin  4>  cos  $  d<p  dA  dt 

The  number  density  in  the  active  zone  of  the  detector 
may  be  found  by  dividing  the  above  expression  by  the  volume 
swept  out  by  the  molecules  in  time  dt . 

_  c^  f  dc  ddi  dA  cos  (j)  d  t 
ns  c  cos  (f)  1  dA'  dt 


(2.2.1) 
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Therefore,  the  number  of  molecules  passing  between 
two  surface  areas,  for  a  gas  at  rest  in  Maxwellian  equili¬ 
brium,  is  found  by  determining  da). 


d  U) 


d  A  ’  cos  (J)  ’ 

9  9 

d  /_  /  c  o  s  z  <J) 


Where  (j>  is  the  angle  between  the  molecular  flight  path  and 
the  normal  of  area  dA .  Substituting  do)  into  expression  2.2.1, 
we  obtain; 


,  c  2  f  d  c  d  A  c  o  s  2  d) 

dnE  -  - ^ - 

since  (p  =  c{)  '  for  this  detector. 

The  numb  er  density  in  the  detector  can  be  d  e t  erm ined 
by  integrating  dns  with  respect  to  c.  Hence, 


2  °°  „  2 
ns  =  ^5-  cos3  *  /  c2  £  dc  =  ^1°  cos3  * 

s  d2  0  4d2 

where  A  =  Tfr  2  .  If  is  small,  then  cos3  41  will  approximately 
equal  unity,  giving 


n  s 


2 

nrr  o 
4d  2 


(2.2.2) 
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In  deriving  expression  2.2.2,  it  was  assumed  that  the 
reflected  velocity  distribution  was  Maxwellian.  Logan  and 
Stickney  (5)  have  shown  that  the  velocity  distribution  varies 
with  the  scattering  angle.  Hinchen  and  Malloy  (6)  compared 
Logan  and  Stickney's  work  to  that  of  Harrison,  who  assumed 
the  reflected  velocity  distribution  was  Maxwellian,  and 
suggested  that  even  if  the  velocity  distribution  is  not 
Maxwellian,  the  error  resulting  from  this  assumption  is  very 
sma 1 1 . 

The  number  of  Molecules  incident  on  the  surface  must 
equal  the  number  of  molecules  reflected.  Hence, 


nicmi 


2  (tt) 


1/2 


nr  cmr 


2  (tt) 


1/2 


where  cm  =  (2RT)1//2  is  the  most  probable  molecular  speed. 

This  reduced  to 


il  =  (ii)2  .  (2.2.3) 

Ti  nr 


Since  the  detector  is  located  in  chamber  B,  there  will 
be  a  background  contribution  to  number  density  due  to  the 
pressure  in  this  chamber.  Therefore,  the  total  number  den¬ 
sity,  n,  monitored  by  the  detector  will  be  the  sum  of  the 
signal,  ng ,  and  the  background,  nbg* 
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n  =  ns  +  nbg  (2. 2. A) 

Substituting  expression  2.2.2  into  expression  2.2.4,  we 
obtain  nr. 


4d' 


nr  =  ~ y  (n_nbg> 


From  expression  2.2.3  it  can  be  shown  that, 


Tr 

Ti 


=  (-  .r°2ni - )2 


4dZ  (n-nbg) 


Hence,  the  thermal  accommodation  coefficient  may  be  written 
a  s  , 


a 


„  r . — ,.-5  P — ^  i. _ i2 

4d^  (n-nb  ) 


(2.2.5) 


and  the  normal  momentum  accommodation  coefficient  may  be 
written  as, 


G  ’ 


1 


4d2  (n-nb  ) 
Ti 


(2.2.6) 
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2 . 3  Working  Conditions 

The  pressure  in  chamber  A  is  kept  constant  at  1  x  10“^ 
t°rr ,  corresponding  to  3.37  x  1 0  ^  3  molecules  per  litre.  From 
conductance  and  pumping  speed  relationships,  the  pressure  in 
chamber  B  will  be  approximately  1.95  x  10“8  torr,  correspond¬ 
ing  to  6.55  x  1011  molecules  per  litre.  The  radius  if  the 
orifice  is  0.1  cm  and  the  distance  d  between  the  orifice  and 

the  active  zone  of  the  detector  is  1.0  cm.  Therefore,  the 

r  2 

geometric  factor,  G  =  — ,  is  2.5  x  10-3. 

4d2 

It  would  be  useful  to  have  upper  and  lox^rer  limits  of 
number  density  in  the  active  zone  of  the  detector  in  order 
to  determine  the  range  of  pressure,  P^,  required  for  the  cali¬ 
bration  of  the  detector.  If  energy  is  completely  accommodated, 
Tr  =  Tg  and  a  =  1.  Thus, 


n  = 


niro' 

4d2 


qA)1'2 

1  s 


+  n 


kg 


The  maximum  value  of  n  will  occur  when  Tg  =  Tj_.  Therefore, 

=  9  x  10^2  molecules  per  litre.  The  minimum  value  of  n 

III  d 

will  occur  when  Tg  is  at  its  maximum  value  of  1293°K.  There- 

12 

fore,  nmin  =  4.65  x  10  molecules  per  litre. 

The  pressure  range  required  in  chamber  A  may  be  deter¬ 
mined  from  the  following  equations. 


r  2  P  a 
L  o  A 

Id2  kTi 


+ 


nb  g 


n 


’ 
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4d2  kT± 

PA  =  - 2 -  (n_nbg) 

'  r  o 

Substituting  the  values  nmax  and  nm^n  for  n,  the  pressure 
range  is  shown  to  be, 

4.7  x  10  ^  1  x  lO-"^  torr. 

The  minimum  value  of  the  signal  background  ratio  will 
be  taken  as  nmin/nbg  which  has  a  value  of  6.1. 

2 . 4  Analysis  of  Experimental  Accuracy 

An  error  analysis  has  been  made  to  determine  the  varia¬ 
tion  in  the  accommodation  coefficient  when  one  parameter  is 
varied  . 

The  square  of  the  standard  deviation  for  a  functional 
relationship  is  defined  as  (7), 


where 


Su2  =  l  (6uj  )2  /n 

j  =  l 


u 


=  (1H)  6xh  +  (|^)  Syj  + 


3x‘ 


3y 


Since  (111)  (111)  will  be  equally  likely  to  be  positive  as 

3x  3y 

o 

negative  for  large  n,  S^ 


becomes 


» 


' 


♦ 


■ 
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s 


u 


2 


(^2 


(|i)2 

dy 


(2.4.1) 


where 


(6x.i)2 

n 


n 


I 


(xj-x)2 

n 


We  are  assuming  all  errors  are  independent  and  may  be  treated 
as  a  normal  distribution. 

Redefining  a  by  grouping  parameters,  we  can  write: 


a 


1  -  [ 


Gn. 


n- 


nbg 


1  - 


Ti 


(2.4.2) 


Since  the  exact  value  for  n  which  corresponds  to  a 
specific  temperature  ratio  is  not  known,  a  value  for  n  and 
Tg/Ti  will  be  chosen  within  the  expected  range.  Thus  n  is 
chosen  as  6.65  x  10J  molecules  per  litre,  and  Ts/T^  as  2.0. 

Applying  expression  2.4.1  to  equation  2.4.2,  the  square 
of  the  standard  deviation  in  a  is, 


. 
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/dct\  2 

9g 


9n 


(fn“)2  S"2 


(2.4.3) 


+ 


(  9a  ^  2 
3nbg 


Differentiating  expression  2.4.2  and  substituting  the 
appropriate  values,  we  get: 


9a 

9G 


2  Gn . 


(n-nbg)2  (1-^) 


=  1.58  x  10 


9a 

9n^ 


2  G  2n  j 


(n-nbg)  (1~0 


=11.7  x  10 


-16 


9a 

3n 


2G2n12 


<"-nbg)  e-ff) 


6.56  x  10 


-13 


3  a 


OP2  2 

2G  nb 


3nbg 


'x  T 

<“-”bg>  (1-Tf 


6.56  x  10 


-13 


9  a 


G^  n  ^ 


~ — s~)  /  (1— 


s^2  _ 


(n-nbg)  T 


(1 


0.97 


. 
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Squaring  each  of  the  above  terms, 


f  3a.  2 


2.5  x  106 


( 9a  )  2 
9n-^ 


1.37  x  10"30 


(2.4.4) 


/ 3a. 2  _  ,  9a  .2 
"  3”bg 


4.3  x  10  2 


C3^-)2  =  0.9  4 

T 

9— 

Ti 


To  determine,  the  sensitivity  of  a  with  respect  to  the 
parameters,  the  standard  deviation  in  the  measured  values 
will  be  assumed  to  be  one  percent  of  their  mean  values.  The 
results  given  in  equations  2.4.6  below  indicate  the  sensiti¬ 
vity  of  a  to  measuring  accuracy  for  each  parameter.  There¬ 
fore,  expression  2.4.1  yields: 

n  n 

l  ( 6x± ) 2  l  (x±-x)2 

S  2  =  J - =  J -  =  (0.01  x)2 

x 


Hence: 


' 
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SG2  =  (6G)2  =  6.25  x  1CT10 


n 


.  =  (6ni ) 2  =  11.35  x  10 

l  1 
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Sn2  =  ( 6 n ) 2  -  44.2  x  1020 


(2.4.5) 


Snbg2=  (««bg)2  =  «.9  x  1018 

ST  2  =  (6~>2  =  4  x  10~4 

s  i  ± 


Substituting  expressions  2.4.4  and  2.4.5  into  expression 
2.4.3,  the  following  square  of  the  standard  deviations  asso¬ 
ciated  with  G,  n^,  n,  n^g  and  Tg/T^  can  be  found. 


,  9  ot  x  2 


=  1.56  x  10 


x  3  oi  \  2 


1.56  x  10 


(i!^)2  s  2  =  1.90  x  10  3 

9n  n 


(2.4.6) 


(~3a  .)2  s„  2  =  1.84  x  10  5 
3nbg'  nbg 
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=  3.76  x  10  ^ 


By  summing  the  expressions  in  2.4.6,  the  square  of  the  stan¬ 
dard  deviation  in  a  can  be  found  (expression  2.4.3).  This 
value  yields  the  standard  error  in  a,  if  all  the  measured 
values  have  an  error  of  one  percent.  Hence,  standard  error 

=  0.0736  ,  fractional  standard  error  =  —  ^an-dar--  error  _  q.092 

a  5 

and  fractional  most  probable  error  =  (.675)  (0.092)  =  0.062 

(fraction  errors  are  based  on  a  =  0.8). 

Since  errors  in  the  measured  parameters  will,  in  all 
likelihood,  not  be  equal  to  one  percent,  the  square  of  the 
standard  deviation  in  a  can  be  determined  by  summing  the  pro¬ 
duct,  (-^? — )2  S  2,  where  S  is  determined  for  each  measured 
9x±  xi  xi 

parameter . 

From  expressions  2.4.6,  it  can  be  seen  that  the  measure¬ 
ment  of  G,  n^,  and  n  are  most  critical.  A  more  rigorous 
analysis  of  the  measurement  of  these  parameters  is  made  in 
the  following  section. 


2 . 5  Evaluation  of  Measurement  Accuracy 

Instrumentation  should  be  chosen  so  that  the  minimal 
experimental  error  can  be  achieved  when  measuring  a  and  a’.- 


2.5.1  Geometry  Applying  expression  2.4.1  to  G  = 

r  2 

-  °-  - ,  we  can  write  the  square  in  the  standard  deviation  in  G 
4d  2 


as  , 


. 
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(2.5.1) 


0.25  x  10~2  Sr  2  +  0.25  x  10“4  Sd2 


o 


Since  rD  and  d  are  determined  by  a  direct  measurement, 
the  square  of  the  standard  deviation  in  these  measurements 
will  be  very  small  providing  extreme  care  is  taken  (a  large 
number  of  actual  measurements  should  be  made).  Since  rQ  is 
machined  and  rigid,  this  measurement  will  have  little  error. 
We  can  assume. 


n 


10  cm2 


(0.0001  r D ) 2 


Due  to  the  difficulty  in  positioning  the  detector  filament, 
the  distance  d  will  have  to  be  measured  a  large  number  of 
times,  and  we  will  assume. 


n 


i 


l  (drd)2 


(0.0001  d)2 


10~8  cm2 


n 


- 

. 
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Substituting  the  appropriate  values  into 
2.5.1,  the  square  of  the  standard  deviation,  S 
to  be  0.5  x  10  1 2 .  The  resulting  uncertainty 
(l^)2  SG2  =  1.25  x  10~2. 


expression 
,  is  found 
in  a  is 


2*5.2.  Temperature  The  measurement  of  temperature 
is  not  critical  except  when  Tg/T^  approaches  unity.  This 
singularity  should  be  noted  and  if  measurements  are  made  in 
this  temperature  range  it  must  be  considered. 

Following  an  analysis  similar  to  that  in  2.5.1,  let 
T  =  Tg/Tj_,  and  the  expression  2.4.1  becomes, 


O 


«  (^1)2 
°  Tg 


+  (1*-) 


-Ts  2 

(— f) 

T  .  ^ 

1  X 


(2.5.2) 


The  temperature  measurements  are  made  with  a  Leeds 
Nor thrup  Potentiometer.  This  millivolt  readout  has  a  resolu¬ 
tion  of  0.005  millivolts  which  corresponds  to  less  than  2°C. 

Therefore  ST  and  ST  can  be  chosen  as  2°C,  and  T  chosen  as 
1  s  1  i 

2  (i.e.,  Ts  =  580°C  and  Tj_  =  290°C).  Substituting  the  above 

9  _ 4 

values  into  expression  2.5.2,  we  find  is  2.38  x  10 

3  a  2  2  -  4 

The  resulting  uncertainty  in  a  is  (-^0  sp  =  2*2^  x  10 


2.5.3. 


Pressure  Gauges 


■ 

' 


■■ 
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2 • j . 3 • 1 «  Pressure  in  Chamb er  A  From  section 
2.4,  it  is  seen  that  the  measurement  of  n^  is  critical. 

The  instrumentation  chosen  to  monitor  this  pressure  must 
have  a  high  degree  of  stability.  For  this  reason,  a  MKS 
Baratron  Model  90H-3  pressure  transducer  was  chosen  to 
monitor  the  pressure  in  chamber  A.  The  manufacturers  claim 
this  gauge  has  a  long  term  drift  of  less  than  one  percent 
after  warming  up. 

One  additional  advantage  of  this  gauge  is  that  it  is 
a  direct  pressure  measuring  device.  It  is  a  diaphram  type 
pressure  transducer  which  is  calibration  against  a  change 
in  capacitance.  Thus,  there  are  no  hot  filaments  to  cause 
degassing  from  the  surrounding  structure  which  may  supply 
impurity  gases  to  contaminate  the  platinum  test  surface. 

2 . 5 . 3 . 2  .  Background  Pressure  A  Varian  partial 
pressure  gauge  model  974-0036,  which  has  a  long  term  drift 
of  less  than  one  percent,  can  be  used  for  monitoring  the 
pressure  in  chamber  B.  The  amounts  of  impurity  gases  pre¬ 
sent  will  also  be  measured  with  this  gauge.  The  length  of 
time  the  test  surface  will  remain  atomically  clean  can  be 
determined  from  the  knowledge  of  the  amounts  of  impurity  gases. 

2.5.4.  Ionization  detector  The  detection  system 
described  in  section  2.1  is  used  for  the  measurement  of  n. 

The  instrumentation  required  for  measuring  the  electron  film 
current  must  have  very  little  long  term  drift.  A  Varian 
dual  range  ionization  gauge  control  unit  model  971-0014  can 


, ; 


' 
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be  adapted  for  this  purpose. 

The  control  unit  grid  is  attached  to  the  detector  shield 
which  corresponds  to  the  electron  collector  (anode) .  The 
potential  of  the  anode  is  +175  volts.  The  filament  shield 
assembly,  which  is  the  electron  emission  side,  is  at  a 
potential  of  +45  volts  and  connected  directly  to  the  control 
unit.  The  ion  collector  is  also  connected  directly  to  the 
control  unit  at  ground  potential.  The  ionization  chamber 
is  set  at  +125  volts  by  a  highly  regulated  Harrison  Labora¬ 
tories  6209B  DC  power  supply. 

A  simple  schematic  for  this  circuit  is  shown  in  Figure 
2.1.1.  The  accelerating  potentials  are  shown  to  provide 
insight  in  the  functioning  of  the  detector.  This  control 
unit  provides  feedback  on  emission  current  which  stabilizes 
it.  This  feedback  is  needed  to  stabilize  the  emmission  cur¬ 
rent  for  accurate  calibration. 

Referring  to  equation  2.2.4  for  the  number  density 
measured  by  the  detector,  we  obtain: 


n  =  Gn±  +  nbg 


(2.5.3) 


and  applying  expression  2.4.1  to  the  above  equation;  the 
actual  square  of  the  standard  deviation  is  obtained, 


2  -  (|i)2  SG2  +  (|T-)2  Sny  +  S, 


n  3  G' 


3nn*  "  ~ni  '3nbg'  nbg 


2  c  _ 2  i  n2 


=  n-j/  SG‘  +  GA  S„  *  +  S 


ni  nb  g 


. 


■ 
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Substituting  the  appropriate  values  from  section  2.4  and 

2.5,  we  find  Sn  =  71.5  x  10^.  Therefore  the  square  of 

the  standard  deviation,  (—^)  ^  S  associated  with  the  un- 

dn  11 

certainty  in  n  is  equal  to  3.07  x  10  ^ .  This  error  is  s  - 
lightly  more  than  an  error  one  percent,  and  it  is  due  to 
measurement  of  n^^  -  in  particular  Sn±  is  too  large. 

Since  n  is  measured  in  terms  of  an  output  current,  the 
error  involved  during  calibration  must  be  considered.  From 
Hagena  (4),  the  equation  for  the  output  current  is  given  as: 

I  =  k  .  ieff  •  leff  •  °eff  •  n  (2.5.4) 

where 

k  is  the  amplification 

ief£  is  the  effective  electron  current 

leff  is  the  effective  electron  flight  path 

aeff  is  the  effective  collision  cross  section 

for  a  nitrogen  molecule  and  an  electron. 

Once  the  system  is  set  in  the  working  condition  described 
previously,  equation  2.5.4  reduces  to: 

I  =  Kn  (2.5.5) 

where  K  is  a  calibration  constant. 

Since  this  calibration  has  not  been  done,  we  will 
assume  a  value  for  K  to  clarify  the  illustration.  We  expect 
the  input  current  to  be  small,  in  the  order  of  nanoamps,  but 
the  degree  of  amplification  will  increase  the  input  current 


. 


.  .  t  \v-  -  [  | :  o 


' 


to  the  order  of  microamps .  Hence,  I  may  be  speculated  to 
be  5  x  10~6  amps  when  n  =  6.65  x  1012  molecules  per  litre. 
Therefore,  k  is  equal  to  7.5  x  10" ^ ^  amp  —  litre. 
Rearranging  expression  2.5.5,  we  have: 


(2.5.6) 


Applying  expression  2.4.1  to  equation  2.5.6 


’K 


2  =  (~)2 


91 


+ 


(!f)2 


n 


n 


W  - 1 


n 


—  9  ■ 

4  n 


where  is  one  percent  of  I.  Substituting  the  correct  val- 
ues,  we  find  the  standard  deviation  in  K  is  1.22  x  10  ,  and 

the  percent  standard  error  is  1.65%.  If  the  points  plotted 
on  the  I  verses  n  graph  are  all  within  +  1.65%  of  the  best 
fit  curve,  then  the  value  of  n  as  determined  during  the  test 
will  have  the  desired  accuracy.  If  the  spread  of  points  on 
the  I  verses  n  plot  is  greater  than  +  1.65%  of  the  best  fit 
curve,  then  adjustments  in  the  control  unit,  and  readout  de¬ 
vice  must  be  made. 

The  critical  measurements  in  the  system  are  n^  and  n. 
The  pressure  in  chamber  A  must  remain  constant  throughout 
the  test  run.  Accurate  pressure  measurements  can  be  obtained 
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with  the  aid  of  a  recorder  which  will  provide  greater  accur¬ 
acy  than  a  panel  meter,  even  if  there  are  no  fluctuations 
in  the  signal  to  be  averaged.  Increased  sensitivity  can 
be  obtained  by  use  of  an  offset  bias  voltage  with  the 
recorder  . 


2 . 6  Experimental  Procedure 

Although  the  number  densities  have  to  be  measured 
accurately  to  obtain  an  accurate  value  of  the  accommodation 
coefficient,  it  will  not  be  necessary  to  have  the  absolute 
value  of  n^  to  the  same  degree  of  accuracy  because  the 

number  density  appears  in  the  equation  for  a  as  a  ratio 

n  • 

- — - .  The  units  of  the  calibration  are  not  important 

n  ~  nb  g 

and  hence  the  absolute  calibration  is  not  critical.  Absolute 

values  of  n^  are  primarily  required  to  determine  if  the 

experiment  is  in  the  free  molecule  flow  regime. 

The  ratio,  _ — _ ,  yields  a  simple  method  for  calibra- 

n  -  nbg 

ting  the  detector.  The  pressure  in  the  chamber  A  may  be 
assumed  correct,  enabling  n  and  nbg  to  be  calibrated  against 
nb  by  using  molecular  flow  theory. 


PB 


C 

c  +  s 


PA 


(2.6.1) 


n 


r02  PA 


4d2  kTi 


+ 


C  +  S  kT 


kg 


(2.6.2) 


may  be  increased  in  increments  from  4  x  10  5  to 


I  • 
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1  x  10  4  torr.  At  each  increment  PA  ,  PB  and  I  should  be 
allowed  to  reach  equilibrium  before  the  readings  are  recorded 
During  calibration  Ts  is  kept  constant  at  the  wall  tempera- 

tur  e ,  . 

A  plot  of  n^  versus  n^gCan  be  made  using  expression 
2.6.1.  Also  from  expression  2.5.5  and  2.6.2,  a  plot  of  I 
versus  n  =  n(PA)  can  be  made  to  determine  K  as  discussed  in 
section  2.5. 

This  calibration  procedure  will  yield  a  method  for 
determining  a,  since  once  I  is  recorded  the  ratio  is  known. 

A  plot  of  I  versus  n-^/n-n^g  can  be  made  for  added  convenience 
In  this  case,  n^  would  be  the  working  number  density  which 
remains  constant  during  the  test  run,  n  and  n^g  are  the 
experimental  values. 

During  the  test  run  the  pressure  in  chamber  A  is  main¬ 
tained  at  1  x  10 “4  torr  with  the  aid  of  a  Granville-Phillips 
automatic  pressure  controller.  The  surface  temperature  is 
increased  and  the  readout  current,  I,  and  pressure,  PB,  are 
recorded  for  each  temperature  station. 

A  plot  of  I  versus  Tg/T-^  can  be  presented  to  show  the 
variation  of  number  density  as  a  function  of  surface  tempera¬ 
ture.  A  plot  of  a  and  a’  versus  Ts/T±  will  exhibit  the  chang 
in  a  and  o'  with  variation  of  the  surface  temperature. 
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CHAPTER  III 


ELECTRON  FLOURESCENCE  DETECTOR 

3 . 1  Description. 

The  electron  beam  technique  used  to  measure  the  number 

density  of  nitrogen,  at  point  P  in  chamber  B  (see  Figure 

¥ 

3.1.2),  is  an  extension  of  a  technique  developed  by  E.O. 
Gadamer  (8).  E.P.  Muntz  (9),  D.J.  Marsden  (10)  and  others 
at  UT IAS . 

When  a  high  energy  electron  collides  with  a  nitrogen 
molecule  which  is  in  the  N2  X^  £g+  ground  state,  the  molecule 
is  ionized  and  electronically  excited  to  the  N2  B  Ju+  state. 
A  photon  of  light  is  emitted  when  the  excited  molecule  makes 
a  spontaneous  transition  to  the  N2+  state.  This 

series  of  transitions  gives  rise  to  the  first  negative  sys¬ 
tem  in  the  nitrogen  emission  spectrum  (9). 

The  electron  gun  used  for  this  detector  is  Phillips 
902-564  television  gun.  The  aiming  system  for  the  electron 
gun  is  designed  to  allow  rotational  and  translational  align¬ 
ments  . 

The  electron  gun  incorporates  electrostatic  lenses  for 
focussing  the  beam.  A  1.5  millimeter  diameter  aperture, 
located  about  3  inches  from  point  P,  is  used  to  collimate 
the  electron  beam. 

The  power  supply  for  operating  the  electron  gun  is  a 
N.J.E.  0-30  KV,  0.0  -  2.0  ma .  The  normal  operating  voltage 


. 
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is  10  KV  ,  and  the  maximum  gross  beam  current  is  about  one 
ma .  The  actual  beam  current  is  monitored  by  a  millimeter. 

The  f lourescence  light  output  is  monitored  by  a  E.M.I. 
9205  S  Photomultiplier.  The  photomultiplier  power  supply 
is  a  Fluke,  model  410B,  0-2.1  KV  which  has  a  high  degree  of 

regulation.  The  photocurrent  is  measured  by  a  Kiethly 
Instruments,  610  B  elec tromenter . 

A  more  detailed  description  of  the  electron  f loure¬ 
scence  detector,  and  its  operation  is  given  in  reference 
(3).  Figure  2.1.1  also  shows  the  basic  optical  features  of 
this  device. 

3 . 2  Theory 

The  light  reaching  the  photomultiplier  is  limited,  to 
the  flourescence  from  the  N2+  (0,0)  band  of  the  nitrogen  ion, 
by  means  of  an  interference  filter  which  only  passes  light 
within  a  narrow  range  of  wavelengths  centred  around  3900 
Angstrom  units. 

The  number  density  of  nitrogen  is  made  up  of  two  parts: 
signal  and  background.  Hence,  we  may  write  the  number  density 
of  nitrogen  at  point  P  (see  Figure  3.2.1)  as, 

n=ns+n^g  (3.2.1) 

where  n„  is  the  contribution  to  the  number  density  of  nitro- 

o 

gen  at  point  P  due  to  molecules  effusing  through  the  orifice 
in  solid  angle  du}g  and  nbg  -*-s  the  contribution  to  the  number 
density  of  nitrogen  at  point  P  due  to  the  partial  pressure 


. 
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of  nitrogen  in  chamber  B.  If  the  number  density  of  the  gas 
in  chamber  B  is  n^ ,  then  n^g  may  be  determined  from  the  mass 
spectrometer  scan  and  the  total  pressure  in  chamber  B.  There¬ 
fore, 


peak  height  of  No+ 

nbe  =  v - —  x  nR  (1-G„) 

®  l  all  peak  heights 

(3.2.2) 

nbg  ~  ^^2  nB 


where  Gg  is  a  geometric  factor. 

The  signal  number  density  is  determined  from  the  follow¬ 
ing  analysis . 


ns 


nn 


db)s 

4tt 


(3.2.3) 


The  solid  angle  dOJ  is  calculated  as  follows: 
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2ir 


d“s  = 


I  / 
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r  dr  d9 
r2  +  d2 


dw  s 


tt  In  ( 
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where  rQ  -  0.1  cm  is  the  radius  of  the  orifice  and  d  =  1  cm 
is  the  perpendicular  distance  between  the  orifice  and  the 
electron  beam.  Substituting  the  appropriate  values,  Gs  is 
equal  to  0.00249.  Hence,  1  -  Gg  is  approximately  equal  to  1.0. 

The  number  density  of  nitrogen  at  point  P  may  be  deter¬ 
mined  by  substituting  expressions  3.2.2,  3.2.3,  and  3.2.4 

into  expression  3,2.1. 

n  ~  nr  FN2  nB  (3.2.5) 

nr  =  q~  (n  fN2  nB ) 

Using  expression  2.2.3,  we  get 


n  -j  o 

(— 

nv 


2 


(3.2.6) 


Substituting  expression  3.2.6  into  expression  1.2.3,  the 
thermal  accommodation  coefficient  may  be  written  as, 


a 


(3.2.7) 


And  from  expression  1.2.4,  the  normal  momentum  accommodation 
coefficient  may  be  written  as, 
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a ' 


1  -  ( 


G  s  n  i 


n  -  FN  n 
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1  -  (~) 
1  i 


1/2 


(3.2.8) 


3 . 3  Working  Condition 

The  upper  and  lower  limits  of  the  number  density  of 
nitrogen  at  point  P  will  enable  us  to  determine  the  required 
calibration  pressure  range  in  chamber  A.  If  energy  is  com¬ 
pletely  accommodated  Tr  =  Tg  and  a  equals  1.0.  From  expres¬ 
sion  3.2.6,  we  can  write: 


n 


Gs  ni(-i)l/2  +  F 


N' 


n 


B 


(3.3.1) 


where  Gg  is  0.00249,  n^  is  3.37  x  10-^  molecules  per  litre, 
and  n^g  =  Fjj  nB  Is  6.55  x  10^  molecules  per  litre.  The 
value  of  n^g  is  only  approximate,  but  can  be  determined 
correctly  when  the  apparatus  is  operating. 

Substituting  the  appropriate  values  into  expression 

-I  O 

3.3.1,  we  find  nmax  =  7.96  x  10  molecules  per  litre  when 
Tr  =  293°K5  and  nm£n  =  4.13  x  10-2  molecules  per  litre  when 
Tr  =  12  9  3  °K . 

Therefore,  the  approximate  pressure  range  required  for- 
calibration  is  , 

-  5  -  4 

4. 7  x  10  <  PA  <  1  x  10  torr 
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The  theoretical  signal  background  ratio  will  be  taken 
as  nmin/nng  and  has  a  value  of  5.3.  However,  this  value 
neglects  stray  light  in  the  system  and  the  dark  current  of 
the  photomult iplier  which  will  make  up  part  of  the  background 
signal . 


3 . 4  Analysis  of  Experimental  Accuracy 

An  error  analysis  similar  to  the  one  completed  in  sec¬ 
tion  2.4  will  be  followed  in  this  section.  Hence,  expression 
2.4.1  may  be  written  as, 


’a 


=  ( 


9a 
9  G, 


+ 

9n-^ 


n. 
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n 


(3.4.1) 


+  (^-)2 

9nbg 


where  n  is  6.65  x  1 0  ^  molecules  per  litre  and  Tg/Tp  is  2.0. 
Values  for  n  and  Tg/Ti  have  been  chosen  within  the  expected 
range . 

Differentiating  expression  3.2.7  with  respect  to  each 
parameter,  we  obtain, 
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1.58  x  103 
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=  11.7  x  10 
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6.56  x  10 
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(n-nbg)3  (I-—-) 


=  6.56x10 


-13 
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G  c  n  0 
1  -  (-2-^-)  2 


n-n 


bg 


T 


s\2 


0.97 


Squaring  each  of  the  above  terms, 


(|“_)2  =2.5  x  106 

3GS' 


c|^-)2 

8n^ 


1.37  x  10  30 


(|^)2 

8n 


(_1°L_)2  =  4.3  x  10  25 
3nbg' 


(3.4.2) 
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To  determine  the  sensitivity  of  a  to  each  parameter, 
the  standard  deviation  in  the  measured  values  will  be  set 
at  one  percent.  Thus,  expression  2.4.1  yields: 
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x 


2 


l 

i  

n 


I  (xi-x)2 

- - - — - —  =  (0.01  x)  2 

n 


hence , 


Sr  2  =  6.25  x  10“10 

Gs 


S2  =  11.35  x  1026 
11  i 


Sn2  =  44.2  x  1020 


Snbg2  =  42.9  x  1018 


Srr.  2  =  4.0  X  10  ^ 

1  S 


T  . 
1 


The  sensitivity  of  0i ,  due  to  variations  in  parameters 
is  determined  by  expression  3.4.1; 


,3a  \  2 
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1.56  x  10 
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3a  v 
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=  1.56  x  10  3 


(f^)2 


n 


1.9  x  10“3 


(3.4.3) 


_9c i_x  2 
3nb  ; 


1.84  x  10~5 


=  3.76  x  10“4 


By  summing  the  expressions  in  3.4.3,  the  square  of  the  stan¬ 
dard  deviation  in  a  can  be  determined  (expression  3.4.1). 

This  value  yields  the  standard  error  in  a,  if  all  the  mea¬ 
sured  values  have  an  error  of  one  percent.  Hence, 

Standard  error  =  0.0552 

Fractional  standard  =  standard_error  =  -069 

a 

error 

Fractional  most  pro-  =  (  0 . 6 75 )  (  0 . 0 6 9 )  =  0.046 

bable  error 

From  expression  3.4.3,  it  can  be  seen  that  the  measure¬ 
ment  of  G,  np  and  n  are  most  critical.  A  more  rigorous  analy¬ 
sis  of  these  parameters  is  made  in  the  following  section. 


3.5  Evaluation  of  Measurement  Accuracy 


Since  the  vacuum  system  used  is  identical  to  the  system 


’ 


36 


used  in  Chapter  II,  some  of  the  instrumentation  associated 
with  the  electron  flourescence  detector  is  identical  to  the 
instrumentation  used  for  the  ionization  detector.  Errors 
associated  with  this  previously  considered  instrumentation 
will  be  repeated  for  completeness. 


3.5.1  Geome try  To  insure  the  correct  error  in  a  is 


known,  a  more  rigorous  treatment  of  the  geometric  factor, 
Gg,  should  be  made.  Applying  expression  2.4.1  to  Gs  =  In 


r  2  +  d2 

( — 2 - ),  we  may  write: 
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)2  s. 
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2d  (r02+d2 ) 


2  sd2 


=  25  x  10“4  Sro2  +  25  x  10  6  Sd2  (3.5.1) 

As  explained  in  section  2.4,  S  r  ^  2  has  a  value  of  appro¬ 
ximately  10~-l'G> 

The  value  of  d  is  not  so  easily  determined  as  in  section 
2.4.  The  distance  d  is  the  average  distance  between  the  ori¬ 
fice  and  the  centre  of  the  electron  beam. 

This  distance  can  be  measured  with  the  aid  of  a  CENCO 
No.  74002  measuring  microscope  made  to  replace  the  micro¬ 
scope  objective  of  a  No.  72905  measuring  microscope. 


This 
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adaption  will  convert  the  microscope  to  a  cat heto meter  with 
a  focusing  range  from  60  cm.  to  infinity.  The  catheto- 
meter  can  be  mounted  on  a  height  gauge  equipped  with  a 
vernier  capable  of  reading  to  0.01  mm. 

The  cathetometer  can  be  focused  at  the  elevation  of 
the  orifice  (bottom  of  chamber  A).  Then  it  can  be  lowered, 
and  focused  on  the  electron  beam.  The  difference  in  the  two 
readings  is  the  distance  d. 

The  pressure  in  chamber  B  will  have  to  be  raised  to 
approximately  10  “  4  torr  before  the  light  intensity  emitted 
from  the  electron  beam  becomes  visible. 

The  centre  of  the  beam  will  be  easily  determined  since 
the  beam  diameter  is  only  1  mm  and  light  is  most  intense  at 
the  centre.  An  estimate  of  the  square  of  the  standard  de¬ 
viation  in  d  is  assumed  to  be  (.001  d)^  cm^  . 

Substituting  values  into  expression  3.3.4  we  find 

Se  2  equals  25.25  x  10  Therefore,  (-^— ) 2  Sg  2  =  6.3  x 

s  d  bg  s 

10“ ^  which  is  the  sensitivity  of  a  due  to  measurement  uncer¬ 
tainties  in  geometry.  This  is  less  than  half  the  value  of 
15.6  x  10“5  obtained  assuming  one  percent  accuracy  in  measure¬ 
ment  s  . 

3.5.2  Temperature  As  described  in  sub-section  2.5.2 

the  square  of  the  standard  deviation  in  c i  associated  with 

X  /X .  is  ( 2  sT  2  =  2.24  x  10“4.  This  error  is  due  to 
s  1  T„  1  s 

3— —  x  . 

T±  1 1 

the  resolution  of  the  Leeds— Nor thrup  Potentiometer. 


. 
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3.5.3  Pressure  Gauges 

3 . 5 . 3 . 1  Pressure  in  Chamber  A  The  square  of 

the  standard  deviation  in  association  with  n,>  is  .)  2  S„  .  2 

1  dn±  ni 

-  1.56  x  1 0 ~ ^ .  This  error  is  due  to  the  one  percent  long 
term  drift  of  the  MK S  Baratron  gauge. 


3.5. 3.2  Background  Pressure  Since  n^g  is 
partial  pressure  of  nitrogen  in  chamber  B,  an  analysis 
be  completed  to  determine  the  error  involved.  Applying 
expression  2.4.1  to  n^g  =  Fjj  nB>  we  can  write> 
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A  value  for  Fj^  (fraction  of  nitrogen  molecules  pre¬ 
sent  in  chamber  B)  may  be  obtained  from  Figure  7  in  reference 
(3).  This  value  is  found  to  be  0.965.  An  estimate  of  the 
standard  error  in  F^^  was  found  from  the  same  plot  by  noting 
the  point  scatter  around  the  best  fit  curve.  This  error  was 


calculated  to  be  1.5  percent  of  the  slope.  Therefore,  S-p 


N 


=  (0.015  FN2)2  and  F^  =  0.965.  Values  for  nB  and  are 

taken  from  section  3.4.  Substituting  into  expression  3.5.2, 

3a  n2  r  2 

nbg 


we 


can  calculate  Sn  2  to  equal  9.2  x  10  and  (^  ^ 
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is  3.96  x  10  2.  This  error  is  not  significant  when  compared 
to  errors  associated  with  the  measurement  of  n^. 

3.5,4  Electron  Flourescence  Detector  This  detection 
system,  which  is  described  in  section  3.1,  is  used  to  moni¬ 
tor  the  number  density  of  nitrogen  at  point  P  ^  description 
of  the  electronics,  operation,  and  calibration  of  this  de¬ 
tection  system  is  given  in  reference  (3). 

Following  a  parallel  analysis  presented  in  sub-section 
2.5.4,  we  can  write  the  equation  for  the  square  of  the  stan¬ 
dard  deviation  in  n.  Therefore,  using  expression  3.2.5, 
we  get : 
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3n 


3n 


bg 


(3.5.3) 


differentiating. 


+  G, 


2 


substituting  the  appropriate  values  from  section  3.4  and  3.5, 
we  find  Sn^  =  7.4  x  1021  and  (-|^-) 2  Sn2  =  3.18  x  10  3.  The 
sensitivity  of  a,  associated  with  the  measurement  of  n  is 
significant  when  compared  to  the  magnitudes  of  the  sensitivi¬ 
ties  for  the  other  parameters.  Sn2  is  large  solely  because 
Sn . 2  is  large,  even  though  Sn^2  can  be  measured  within  a  one 
percent  standard  error. 
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Since  n  is  measured  in  terms  of  an  output  current, 
the  error  involved  during  calibration  must  be  considered. 

In  reference  (10),  the  emission  intensity  is  given  by, 

I'  =  n’  he  V  Anm  (3.5.4) 

where 

n'  =  the  number  density  of  the  molecules  in  the 
excited  state  at  any  instant 

h  =  planck's  constant 

c  =  velocity  of  light 

V  =  the  wave  number  of  the  light  emitted 

Anm  =  Einstein's  spontaneous  emission  probability 

The  depopulation  of  the  excited  state  by  spontaneously 
emitting  radiation  per  second  is  n'Anm  and  the  depopulation 
of  the  excited  state  by  the  collisions  with  neutral  mole¬ 
cules  per  second  is  2n'n  62  (4ttRT)^^  (see  reference  (10)). 

Hence , 

F.n  =  n'Anm  +  2n'n  62  (4TTRT)1/2  (3.5.5) 

where , 

n  =  number  density  of  the  gas 

F  =  fraction  of  n  excited  by  the  electron  beam 
per  second 

6  = 

There  fore,. 


effective  radius  of  the  collision  cross  section 


. 
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n'  = 


F  .n 


Anm  +  2n  62  1/2 


(3.5.6) 


Substituting  equation  3.5.6  into  equation  3.5.4,  we  can  write 


I'  = 


F  .  n  he 


(1  +  2 n  6 2  (4nRT) 1/2/Anm) 


(3.5.7) 


The  quenching  collision  term  is  not  significant  if 


2n  6 2  (4ttRT)  1//2/Anm  <<  1 


and  since  n  =  6.65  x  1012  per  litre,  the  quenching  term  is 
not  significant.  Therefore, 


or 


I '  =  F  .  h  .  c  .  n 

I  =  k  n 


(3.5.8) 


where , 

I  =  the  output  current  of  the  photomultiplier 
K  =  a  calibration  constant 

This  constant  has  been  determine  in  reference  (3),  which  also 
cites  the  factors  incorporated  in  it.  The  value  of  K  was 
found  to  be  4.1  x  10“5  Amp/torr  which  corresponds  to  1.24  x 
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10  2^  Amp-litre. 

Folio w in g  a  parallel  analysis  as  presented  in  sub-sec¬ 
tion  2.4.3,  the  desired  accuracy  in  K  can  be  found.  Apply- 
ing  expression  2.4.1  to  K  =  1/n,  we  can  write, 
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(M)2 

V3n 


(3K)2 

dl' 


(3.5.9) 


where , 

n  =  6.65  x  10^/litre 

I  =  8.2  x  10~12  ^mp 

Sj2  =  (.01  I)2  =  67.2  x  lO"28  Amp2 


we  have  Sjr2  =  4  x  10  -*2 


»  ~K 


Substituting  the  appropriate  values  into  expression  3.5.9, 

Sv  =  2  x  10-2^  Amp-litre,  and  the 
percent  standard  error  is  1.6%.  Therefore,  n  can  be  deter¬ 
mined  accurately,  if  the  spread  of  points  on  the  I  versus  n 
plot  is  not  greater  than  +  1.6  percent  of  the  best  fit  curve. 
The  critical  measurement  for  this  system  is  n-^.  The 

smaller  the  standard  error  in  n^,  the  smaller  the  total  error 

3  cx  2  o 

in  a.  The  largest  contributors  to  the  error  in  a  is  (-— — )  Sn  . 

dnj_  “x 

Pi  2  2 

and  (iL9b  S„  (S„  is  dependent  upon  n  •  )  .  Since  the  measurement 
3n  n  n  1 

of  n^  is  the  most  important  factor  in  determining  the  accuracy 
of  a,  a  recorder  should  be  used  as  the  output  instrument  n^ 
to  detect  drift  and  average  the  signal  over  a  period  of  time. 


. 
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It  will  enable  the  investigator  to  determine  Sn  more  pre- 

11  i 

cisely . 

3 . 6  Experimental  Procedure 

The  electron  flourescence  device  can  be  calibrated 
with  respect  to  n^  .  Because  of  the  number  density  ratio, 
the  absolute  value  of  n  ^  does  not  have  to  be  known  as  accur¬ 
ately  as  the  ratio. 

The  permanent  magnet  on  the  partial  pressure  gauge 
deflects  the  electron  beam,  and  has  to  be  removed  prior  to 
the  tests.  The  partial  pressure  of  nitrogen  is  obtained  from 
a  calibration  of  partial  pressure  of  nitrogen  against  total 
pressure,  made  prior  to  operating  the  electron  beam. 

After  the  above  calibration  has  been  completed,  the 
test  surface  can  be  cleaned,  and  the  water  circulated  through 
the  cooling  jacket  in  chamber  A. 

The  output  current,  I,  from  the  photomultiplier  is  then 
calibrated  against  n^,  and  by  using  expression  3.2.5  it  can 
be  calibrated  against  n.  From  the  plot  of  I  versus  n,  the 
value  of  K  and  can  be  determined.  The  percent  standard 

error  in  K  must  not  be  greater  than  1.6%,  if  the  desired 
accuracy  for  the  measurement  of  n  is  to  be  attained. 

A  plot  of  I  versus  n^/n-n^g  can  be  made.  In  this  case, 
n is  the  operating  number  density  which  remains  constant 
throughout  the  experiment;  n  and  n^g  are  the  experimental 


values . 
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Once  the  system  has  been  calibrated  no  adjustment  of 
the  instrumentation  can  be  tolerated. 


CHAPTER  IV 


DISCUSSION 

4.1  Comparison  of  Detectors 

In  both  the  ionization  and  electron  f lourescence  detec¬ 
tors,  the  major  source  of  error  occurs  in  the  measurement  of 
n^ .  It  is  imperative  that  Sn <  be  minimized.  A  recorder  can 

be  used  to  average  the  output.  This  will  reduce  S  provided 

11  i 

there  is  no  long  term  drift,  either  in  the  pressure  or  in 
the  instrumentation. 

The  advantages  of  the  ionization  detector  are;  direct 
measurement  of  number  density,  and  signal  background  ratio 
(nmin/nbg)  i  s  aPProxiniately  six. 

The  limiting  factor  is  the  x-ray  limit  of  the  ioniza¬ 
tion  detector,  and  the  background  number  density  in  chamber  B. 

The  advantages  of  the  electron  flourescence  detector 
are:  the  number  density  of  nitrogen  can  be  measured  indepen¬ 
dently  of  impurity  gases,  and  the  electron  beam  does  not  dis¬ 
turb  the  molecular  flight  path. 

The  disadvantages  of  this  device  are:  the  signal  inten¬ 
sity  is  low  because  photons  emitted  are  collected  by  the 
photomultiplier  and  only  a  limited  number  of  gases  have 
emission  spectra  of  sufficient  intensity  to  be  useable. 

4 . 2  Modifications  to  the  System 

The  baffle  obstructing  the  conductance  between  chambers 


' 
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B  and  C  should  be  removed  to  insure  that  the  effective  pump¬ 
ing  speed  from  chamber  B  is  1900  litres  per  second.  Two 
smaller  baffles  can  be  placed  over  the  titanium  filament 
and  the  ion  pump. 

The  MKS  Baritron  gauge  has  proven  to  be  troublesome. 
Since  the  measurement  of  n^  is  most  important,  it  should  be 
replaced  by  a  more  reliable  gauge.  A  Varian  Dual  Range 
Ionization  gauge  wh ich  has  a  thoria-coated  iridium  filament 
(low  work  function),  could  be  used  to  replace  the  Baritron 
gauge.  This  control  unit  has  a  long  term  drift  of  less  than 
one  percent . 

The  ionization  detector  is  functional  at  the  present 
time.  However,  a  slight  modification  on  the  filament  mounts 
would  improve  the  convenience  of  replacing  them. 

The  electron  beam  system  should  have  the  change  in 
pressure  gauge  cited  above  as  well  as  the  following  changes: 

1.  Light  baffles  could  be  placed  around  the  electron 
beam  except  at  point  P.  Baffles  can  reduce  the  photocurrent 
due  to  the  excitation  of  the  background  molecules. 

2.  The  orifice  could  be  moved  closer  to  the  electron 
beam.  If  the  distance  d  is  reduced  by  half,  the  signal  num¬ 
ber  density  will  subsequently  be  increased  by  four,  and  the 
photocurrent  by  the  same  amount. 

3.  If  the  effective  electron  beam  current  could  be 
increased  two  fold,  the  photomultiplier  output  current  will 
be  increased  proportionally. 

4.  A  concave  mirror  could  be  placed  opposite  the  lens 


. 
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and  point  P.  Point  P  can  be  placed  at  the  focal  point  of 
the  lens  in  the  image  of  the  mirror.  This  mirror  will  in¬ 
crease  the  effective  solid  angle  of  light  received  by  the 
photomultiplier.  For  example  a  mirror  with  an  effective  area 
of  25  square  inches  placed  2  inches  from  point  P  would  result 
in  an  increase  in  solid  angle  of  approximately  24  times  the 
solid  angle  obtained  in  reference  (3) . 

These  modifications  will  increase  the  calibration 
constant  approximately  190  times  the  previous  value.  There¬ 
fore,  the  expected  value  of  K  would  be  2.3  x  10“22  amp-litre, 
and  the  expected  photomultiplier  current  would  be  1.6  x  10-^ 
amps,  which. is  three  orders  of  magnitude  larger  than  the 
photomultiplier  dark  current. 

The  analysis  presented  in  Chapters  II  and  III  is  only 
approximate,  since  the  actual  value'  for  n  corresponding  to 
a  given  temperature  Ts  is  not  known.  The  actual  values  can 
not  be  known  until  the  complete  system  is  calibrated.  There¬ 
fore,  all  the  errors  previously  calculated  are  subjected  to 
the  dependence  of  n  on  Ts.  However,  these  variations  are 
not  likely  to  be  large. 

Manufacturer’s  specifications  for  the  control  units 
previously  mentioned  all  possess  a  long  term  drift  of  less 
than  one  percent.  For  any  given  instrument,  the  error 
associated  with  the  measurement  of  np  will  probably  be  less' 
than  the  error  calculated  in  this  study.  Therefore,  the 
actual  error  in  a  could  be  smaller  than  the  error  estimated 


in  this  study. 
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CHAPTER  V 


CONCLUSIONS 


When  using  the  ionization  detector  to  measure  a,  the 
following  sensitivities  are  associated  with  uncertainties 
in  the  parameters: 

(||)2  SG2  =  0.001  X  10-3 

(l^)2  S_  2  =  1.56  X  10-3 

3n^  11  i 

(3ot)2  s  2  =  3.07  x  10"3 

dn  11 

(  3o^_)2  s  2  =  Q.02  x  10“3 
9*bg  nbg 

(3o_)2  s  2  =  o.22  x  10“3 

T  1  s 


2 

The  square  of  the  standard  deviation,  Sa  ,  is  A. 86  x 
10-3.  The  percent  standard  error  in  a  is  8.7%  and  the  per 

cent  most  probable  error  in  ot  is  5.9%. 

The  square  of  the  standard  deviation  in  o  (as  deter¬ 
mined  from  the  ionization  detector  error  analysis)  is  3.1 
10"3.  The  percent  standard  error  is  5.0%  and  the  percent 


' 
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most  probable  error  is  4.7%. 

The  sensitivity  of  a  has  the  following  dependence  when 
the  electron  flourescence  detector  is  used: 


( 


3a  >. 
3GS; 


2 


/  3a 
dn  ^ 


) 


2 


(|^)2 


=  0.06  x  10  3 


=  1.56  x  10  3 


=  3.18  x  10“3 


^L-)2  S„  2 
8nbg 


0.003  x  10~3 


( 


2 


=0.22  x  10" 3 


The  square  of  the  standard  deviation  in  a  is  5.02  x  10“3 
The  percent  standard  error  in  a  is  8.9%,  and  the  percent 
most  probable  error  is  6.0%. 

The  square  of  the  standard  deviation  in  a',  as  deter¬ 
mined  from  the  electron  flourescence  detector  error  analysis, 
is  3.64  x  10~3.  The  percent  standard  error  is  7.5%  and  the 
percent  most  probable  error  is  5.0%. 

This  type  of  system  is  well  suited  for  the  measurement 
of  thermal  and  normal  momentum  accommodation  coefficients. 
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The  experimental  accuracy  attained  by  using  either  detection 
system  is  approximately  six  percent.  The  major  source  of 
error  occurs  in  the  measurement  of  n^ . 
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SCHEMATIC  CROSS  SECTION  OF  IONIZATION  DETECTOR 
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SCHEMATIC  OF  GEOMETRICAL  ARRANGEMENT  NEAR  TARGET 
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SCHEMATIC  PLAN  VIEW  OF  ELECTRON  BEAM  AND  OPTICAL  SYSTEM 
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